Melting and freezing of argon in a granular packing of linear mesopore arrays 
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Freezing and melting of Ar condensed in a granular packing of template-grown arrays of linear 
mesopores (SBA-15, mean pore diameter 8 nm) has been studied by specific heat measurements C 
as a function of fractional filling of the pores. While interfacial melting leads to a single melting 
peak in C , homogeneous and heterogeneous freezing along with a delayering transition for partial 
fillings of the pores result in a complex freezing mechanism explainable only by a consideration of 
regular adsorption sites (in the cylindrical mesopores) and irregular adsorption sites (in niches of 
the rough external surfaces of the grains, and at points of mutual contact of the powder grains). 
The tensile pressure release upon reaching bulk liquid/vapor coexistence quantitatively accounts for 
an upward shift of the melting/freeezing temperature observed while overfilling the mesopores. 

PACS numbers: 64.70.Nd, 65.80.+n, 65.40.Ba 



The arguably most conspicuous effect of pore confine- 
ment on molecular condensates is the shift of phase tran- 
sitions with respect to the non-confined bulk state. Thus 
condensation (the vapor-liquid transition) does not occur 
at the saturated vapor pressure po of the bulk liquid but 
in a wide range of reduced vapor pressures P = p/po, 
P < 1, starting with the condensation of an adsorbed 
layer on the pore walls and culminating in pore filling via 
capillary condensation[l|. In a similar way melting and 
freezing of the material in the pores occurs at tempera- 
tures T well below the bulk triple point temperature T3 
Q. These reductions are usually interpreted in terms of 
the Kelvin- and the Gibbs-Kelvin equation, resp. Here 
interfacial energies are introduced and because of their 
competition with volume free energies, the shift scales 
with the inverse of the characteristic linear dimension L 
of the geometric configuration. In case of tubular pores 
this is the pore radius R. 

Here we present an experimental calorimctric study on 
the melting and freezing of one of the most simple con- 
densates imaginable, the rare gas Ar, in a template grown 
mcsoporous silica, known as SBA-15 and considered one 
of the most ideal mesoporous substrates available 0]. 

The preparation has been described elsewhere [H • The 
material obtained is a powder of grains, about 1 micron in 
size, which have a relatively rough surface (Fig. Q}. The 
pores within the grains are linear, non-ramified, parallel, 
have a relatively uniform cross section, and form a peri- 
odic 2D hexagonal array. Hence network effects and the 
blocking of the solidification front in bottle necks [j| are 
not expected to play a major role. The Bragg angles of 
the diffraction pattern of the empty pore lattice give a 
lattice parameter (= porc-porc-distance) of 10.7 nm 
A fit of a radial electron density p(r) profile to the Bragg 
intensities (with p = for r < Rq — d , p = constant for 
r > Rq + d, and with a linear function in the " corona" in 
between yields a pore radius Rq of 3.8 nm, correspond- 
ing to a porosity of 0.49, and a corona thickness 2d of 



FIG. 1: Ar sorption isotherm in SBA-15, measured at 86 K. 
Inset: Scanning electron micrographs of a SBA-15 grain taken 
with different spatial resolutions as indicated in the figure. 



2.2 nm. A volumetric Ar sorption isotherm (normal- 
ized uptake / = n/no as function of the reduced vapour 
pressure P, where n is the adsorbed amount of Ar and 
no = 30 mmol is the amount of Ar for complete filling of 
the pores), recorded in the liquid regime of the Ar pore 
filling at 86 K, is shown in Fig. [1] The initial reversible 
part is due to the condensation of Ar onto the pore walls, 
the hysteretic part to the filling of the pores by capillary 
condensation on p-increase (starting at about /a = 0.65) 
and to the evaporation of the capillary condensate on p- 
decrease (terminating at about /d = 0.42). The branches 
of the hysteresis loop arc relatively steep compared to 
what is observed for other mcsoporous substrates, sug- 
gesting a narrow distribution of the pore radius and the 
absence of pore network effects such as pore blocking [(| . 
The final slope in the post-filling regime / > 1 (absent 
e.g. for monolithic Vycor substrates) demonstrates that 
there are further sites available for condensation, with 
characteristic linear dimensions larger than Rq, but still 
finite. We think of tapered pore mouths, niches of the 
rough external surface of the grains, and the points at 
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which the powder grains arc in mutual contact. 
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FIG. 2: (color online). Specific heat curves of cooling (solid 
line) and heating (dashed line) cycles for selected filling frac- 
tions / indicated in the figure. 

The calorimeter consists of the sample cell, heat shield, 
and the outer jacket, with two separate vacuum spaces 
in between. The sample cell and the heat shield are 
equipped with thermometers (1000 fl Pt-resistors) and 
resistance heaters. The heat shield is in thermal contact 
with the cold plate of a closed cycle He refrigerator. 

The cell contains 1.3 g SBA-15 and a coil of Ag wire 
(that is meant to reduce the thermal time constant within 
the powder). 12 fractional fillings have been investigated 
that have been prepared by introducing suitable amounts 
of Ar into the cell via a fill line at 86 K. 

We employ a scanning experiment the calorimetric sig- 
nal is proportional to dQ/dT, the amount of heat per unit 
T-interval that goes into or comes out of the sample. 
dQ/dT = KAT/r, AT and K arc the temperature dif- 
ference and the heat leak between cell and shield, r is 
the cooling/heating rate of the shield. The heat leak has 
been realized by a He gas pressure of a few mbar in the 
vacuum space between cell and shield, AT was typically 
1 K, the heating/cooling rate was held constant at a value 
of 0.05 K/min. Quite elaborate thermometry allowed us 
to work with such low rates which are more than one 
order of magnitude smaller than what is standard with 
commercial equipment. 

The samples have been subject to a complete freez- 
ing/heating cycle starting at 86 K in the liquid regime, 
taking the sample to a minimum temperature of 50 K 



deep in the solid regime, and bringing it back to 86 K. 
The results for a selected set of fillings are shown sepa- 
rately for cooling and heating in Fig. [2] 

Samples with / > 0.45 do show heating and freez- 
ing anomalies, samples with / < 0.38 do not. Analo- 
gous observations have been made for Ar in Vycor jjjj]. 
Thus a phase-transition like change from the liquid to 
the crystalline state is basically reserved to the capillary 
condensed part of the pore filling. The initial "dead" 
fraction of the condensate up to about / = 0.4, adsorbed 
on the pore walls and filling the niches of the corona, does 
not participate in a collective freezing melting process. 

Melting of the "mobile" part occurs in a relatively 
narrow T-interval of 2 K, which for / < 1 is centred 
at 74.5 K, whereas freezing - say for / = 0.92 - ex- 
tends over 7 K from 71.5 K down to 64.5 K with the 
freezing anomaly having a rather complex shape (Fig. 
[2]). The entropy of fusion, So, is obtained by integrat- 
ing the calorimetric signal above the smooth background, 
So = c J{dQ'/dT)/TdT. The coefficient c has not been 
determined, but is the same for all heating and cooling 
runs. A plot of So as function of / is shown in Fig. [3] 
The values of So obtained from cooling and heating runs 
are identical within experimental error. This is not only 
a consistency check but also means that all parts of the 
freezing anomalies are indeed due to no other phase tran- 
sition but solidification. For /< 1, So is a linear function 
of / that extrapolates to zero at about / = 0.39, close to 
/d of the 86 K-sorption isotherm. Thus the part of the 
condensate that supports the collective freezing/melting 
process comprises not only what is formed by capillary 
condensation but also the fraction which grows initially 
in form of a metastabel film at the pore walls (/a — /d)- 

For ease of discussion we divide the freezing anomaly 
into the parts I, II, III which are meant to pertain to cer- 
tain fractions of the pore filling. Part II and III are the 
peaks at 67.5 K and 65.5 K, respectively, part I the re- 
mainder at higher T. The /-dependence of the three 
parts is also shown in Fig. [3l 

The fact that the pore filling solidifies sequentially, but 
melts practically in a single step strongly suggests that 
solidification leads to a reorganization of the capillary 
condensate upon which the different fractions I, II, III of 
pore material loose their identity. 

The melting anomalies have an asymmetric shape. For 
/ < 1 , the low-T wings of the anomalies for different fill- 
ing fractions superimpose. We have interpreted the anal- 
ogous observations on Ar in Vycor in terms of interfacial 
melting in a cylindrical pore starting at the boundary 
between the dead and the mobile part of the pore filling 
rather than by referring to a pore size distribution which 
had to be rather special in order to reproduce the pecu- 
liar shape of the anomalies 0]. This interfacial melting 
scenario is similar to surface melting at semi-confined, 
planar surfaces Q. Calorimetry can of course not prove 
whether the melting front propagates along the pores or 
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radially inward, but the mere fact that the asymmetric 
shape of the anomaly is recovered in the highly homoge- 
neous pores of the present substrate is a strong argument 
in favor of intcrfacial i.e. radial melting. 

As for freezing, part III of the freezing anomaly rep- 
resents the component of the filling that is the last to 
solidify on cooling. So (III) shows the peculiar type of 
/-dependence , So (III) scales with the vapor filled frac- 
tion of pore space. So (III) is due to the solidification of 
the mobile part of the liquid film on the pore walls in 
otherwise empty pore segments. Upon solidification the 
film delayers and the material involved joins the already 
solidified capillary condensate (and melts as such on sub- 
sequent heating) . The same observations have been made 
for Ar in Vycor 0]. 

Still referring to the underfilled situation, / < 1, this 
reasoning leaves parts I and II for the freezing of the 
capillary condensate in the regular pores. We recall that 
the intcrfacial melting model not only explains the low- 
ering of the equilibrium freezing/melting temperature To 
in pore confinement but also predicts metastable states, 
liquid ones down to the lower spinodal temperature T~ 
and solid ones up to the upper spinodal temperature T + . 
We propose to identify peak II with the freezing of the 
capillary condensate via homogeneous nucleation of the 
solid phase at T~ . Peak II is sharp because of the narrow 
pore size distribution of SBA-15. Part I of the freezing 
anomaly on the other hand stems from the freezing via 
heterogeneous nucleation within the rather broad (To, 
T~ )-interval, triggered by the presence of various types 
of small nuclei residing at irregular sites. The chance that 
a parcel of liquid is in contact with such nuclei increases 
with /. That is why part II dominates (in coexistence 
with part III) at low / and goes through a maximum 
slightly below / = 1 whereas part I gradually develops 
with increasing / in a somewhat delayed fashion (Figs. [2] 
and EJ). 

The data on the overfilled sample, 1.04 < / < 1.19, 
shows that the freezing/melting of different parts of the 
Ar condensate are not independent from each other, an 
observation that supports the idea of heterogeneous freez- 
ing. Quite in contrast to the results on a monolithic Vy- 
cor sample, the melting and freezing of the excess mate- 
rial, / — 1, does not lead to a (5-peak at the triple point of 
the bulk system (at 83.8 K), but in first respect to some 
extra calorimetric signal at the high-T end of the anoma- 
lies, both for freezing and melting. Obviously most of the 
excess material resides in sites, say in niches of the exter- 
nal surface of the powder grains, with dimensions that 
are still in the nm-range just slightly larger than the di- 
ameter of the regular pores. The presence of this material 
has an effect on the freezing and melting of the material 
in the pores. The melting peak for / = 0.92 for example 
is centered at T m = 74.8 K, but adding excess material 
the melting peak eventually shifts to 76.1 K and even 
reduces the calorimetric signal at the former T m . Thus 
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FIG. 3: Entropy plots of the different anomaly parts as a 
function of filling fraction /. The dashed lines are guides for 
the eye. 



the presence of material outside stabilizes the solid state 
in the pores. 

The pore filling experiences a negative hydrostatic 
pressure p^ of 70 ± 10 bar for / < 1 Q. This value 
is dictated by the Laplace equation while assuming semi- 
spherical, concave menisci with negative curvature radii 
of the condensate- vapor interface on the order of the pore 
radius [Icj . Upon post-filling of the pores the negative 
curvature of these interfaces gradually vanishes; hence 
the tensile pressure in the pore condensate disappears. 
Based on the p^ dependency of the melting line of Ar [llj , 
an increase of the melting temperature of 1.7 ± 0.3 K is 
expected for a 70±10 pressure release, in good agreement 
with the ~ 1.5K T-shift of the melting peak documented 
in Fig. 2. The upward shift of the onset of the freezing 
anomaly for / > 1 (part I) can be explained in an anal- 
ogous way. An opposite, downward shift of T m upon 
approaching / = 1, observed for water in mesopores [12| . 
is consistent with our tensile pressure hypothesis due to 
water's anomalous T m (ph)- curve. 

For / = 0.66 partial cooling/heating cycles have been 
investigated (Fig. In the first type of such cycles a 
part of the condensate is solidified by cooling to some 
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FIG. 4: (color online). Incomplete (partial) cooling-heating 
cycles at / = 0.66 (solid lines) in comparison to heating, 
cooling curves after complete solidification, melting, resp. at 
/ = 0.62 (dotted lines). 



temperature within the T-range of the freezing (Fig. 0}t) 
and the melting of this part is then studied on subsequent 
heating (Fig. [4Jd). In the first run the freezing of part I 
is initiated by cooling down to 68 K, in the second one 
parts I and II are solidified by cooling down to 66.5 K. 
As far as position and shape are concerned, the resulting 
melting anomalies are identical to what is obtained after 
complete solidification (see Fig. [2] for data on / = 0.62 
and 0.72), it is just that So is scaled down in proportion to 
the amount of material solidified on cooling. Parts of the 
condensate with largely different freezing temperatures 
melt in the same narrow T-interval. 

The second type of cycles starts at 60 K deep in the 
solid regime and a part of the solid is melted by heat- 
ing up to temperatures within the T-range of melting 
(Fig. [4]i) . The re-solidification is then studied by subse- 
quent cooling (Fig. [4jl). In such partial cycles the calori- 
metric signal right at the onset of freezing is enhanced. 
Obviously the fraction of material that has been melted 
in the heating cycle is still in contact with the solid re- 
maining and can therefore resolidify directly without hav- 
ing to overcome a nuclcation barrier. Peak III is recov- 
ered in the partial cycles. Whenever some part of the 
pore filling is liquefied, the mobile part of film coating 
on the pore walls is re-established. As for peak II, the 
completed freezing process of the first cooling run down 
to 60 K leads to a re-arrangement of the pore filling in 
pore space by which the isolated parcels are eliminated, 
the freezing of which (via homogeneous nucleation) gives 
rise to peak II. These parcels are not re-established in 
case melting along the heating run is incomplete. This is 
why peak II is absent in the lower trace of Fig. 4d. Even 
in case of complete melting by heating up to a temper- 
ature of 75.5 K slightly above the high-T cutoff of the 
melting anomaly, the isolated parcels are not fully recov- 



ered compared to the situation of the virgin sample right 
after preparation by vapor condensation at 86 K, see the 
dotted line in Fig. 4. This explains the reduced size of 
peak II in the upper trace of Fig. 4d. After a freezing- 
melting cycle the pore liquid keeps some memory of the 
spatial distribution of the pore solid 
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The freezing/melting phenomenology encountered in 
SBA-15 challenges the common notion which underlies 
pore size spectroscopy techniques ("thcrmoporometry") 
[14j relying on simple geometric relations between pore 
diameter and specific heat anomalies. As demonstrated 
here for a prototype mesoporous system, neither a re- 
sorting to the specific heat anomaly in freezing nor in 
melting will help in a granular packing of mesoporous 
grains. Freezing is demonstrated to be significantly af- 
fected by heterogeneous nucleation, not to mention the 
delayering transition for partial fillings, whereas melting 
is governed by intcrfacial melting processes. 

Freezing and melting in such a granular packing of 
comparably simple linear mcsoporcs is clearly a compli- 
cated process. Metastable states, an interplay of homo- 
geneous and heterogeneous nucleation processes and the 
coupling between different regions of pore space lead to 
a remarkable complex phenomenology. 
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